home *** CD-ROM | disk | FTP | other *** search
/ IRIX Base Documentation 2002 November / SGI IRIX Base Documentation 2002 November.iso / usr / share / catman / p_man / cat3 / SCSL / clatzm.z / clatzm
Encoding:
Text File  |  2002-10-03  |  4.1 KB  |  133 lines
  1.  
  2.  
  3.  
  4. CCCCLLLLAAAATTTTZZZZMMMM((((3333SSSS))))                                                          CCCCLLLLAAAATTTTZZZZMMMM((((3333SSSS))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      CLATZM - routine is deprecated and has been replaced by routine CUNMRZ
  10.  
  11. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  12.      SUBROUTINE CLATZM( SIDE, M, N, V, INCV, TAU, C1, C2, LDC, WORK )
  13.  
  14.          CHARACTER      SIDE
  15.  
  16.          INTEGER        INCV, LDC, M, N
  17.  
  18.          COMPLEX        TAU
  19.  
  20.          COMPLEX        C1( LDC, * ), C2( LDC, * ), V( * ), WORK( * )
  21.  
  22. IIIIMMMMPPPPLLLLEEEEMMMMEEEENNNNTTTTAAAATTTTIIIIOOOONNNN
  23.      These routines are part of the SCSL Scientific Library and can be loaded
  24.      using either the -lscs or the -lscs_mp option.  The -lscs_mp option
  25.      directs the linker to use the multi-processor version of the library.
  26.  
  27.      When linking to SCSL with -lscs or -lscs_mp, the default integer size is
  28.      4 bytes (32 bits). Another version of SCSL is available in which integers
  29.      are 8 bytes (64 bits).  This version allows the user access to larger
  30.      memory sizes and helps when porting legacy Cray codes.  It can be loaded
  31.      by using the -lscs_i8 option or the -lscs_i8_mp option. A program may use
  32.      only one of the two versions; 4-byte integer and 8-byte integer library
  33.      calls cannot be mixed.
  34.  
  35. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  36.      This routine is deprecated and has been replaced by routine CUNMRZ.
  37.      CLATZM applies a Householder matrix generated by CTZRQF to a matrix.
  38.  
  39.      Let P = I - tau*u*u',   u = ( 1 ),
  40.                                  ( v )
  41.      where v is an (m-1) vector if SIDE = 'L', or a (n-1) vector if SIDE =
  42.      'R'.
  43.  
  44.      If SIDE equals 'L', let
  45.             C = [ C1 ] 1
  46.                 [ C2 ] m-1
  47.                   n
  48.      Then C is overwritten by P*C.
  49.  
  50.      If SIDE equals 'R', let
  51.             C = [ C1, C2 ] m
  52.                    1  n-1
  53.      Then C is overwritten by C*P.
  54.  
  55.  
  56. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  57.      SIDE    (input) CHARACTER*1
  58.              = 'L': form P * C
  59.              = 'R': form C * P
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. CCCCLLLLAAAATTTTZZZZMMMM((((3333SSSS))))                                                          CCCCLLLLAAAATTTTZZZZMMMM((((3333SSSS))))
  71.  
  72.  
  73.  
  74.      M       (input) INTEGER
  75.              The number of rows of the matrix C.
  76.  
  77.      N       (input) INTEGER
  78.              The number of columns of the matrix C.
  79.  
  80.      V       (input) COMPLEX array, dimension
  81.              (1 + (M-1)*abs(INCV)) if SIDE = 'L' (1 + (N-1)*abs(INCV)) if SIDE
  82.              = 'R' The vector v in the representation of P. V is not used if
  83.              TAU = 0.
  84.  
  85.      INCV    (input) INTEGER
  86.              The increment between elements of v. INCV <> 0
  87.  
  88.      TAU     (input) COMPLEX
  89.              The value tau in the representation of P.
  90.  
  91.      C1      (input/output) COMPLEX array, dimension
  92.              (LDC,N) if SIDE = 'L' (M,1)   if SIDE = 'R' On entry, the n-
  93.              vector C1 if SIDE = 'L', or the m-vector C1 if SIDE = 'R'.
  94.  
  95.              On exit, the first row of P*C if SIDE = 'L', or the first column
  96.              of C*P if SIDE = 'R'.
  97.  
  98.      C2      (input/output) COMPLEX array, dimension
  99.              (LDC, N)   if SIDE = 'L' (LDC, N-1) if SIDE = 'R' On entry, the
  100.              (m - 1) x n matrix C2 if SIDE = 'L', or the m x (n - 1) matrix C2
  101.              if SIDE = 'R'.
  102.  
  103.              On exit, rows 2:m of P*C if SIDE = 'L', or columns 2:m of C*P if
  104.              SIDE = 'R'.
  105.  
  106.      LDC     (input) INTEGER
  107.              The leading dimension of the arrays C1 and C2.  LDC >= max(1,M).
  108.  
  109.      WORK    (workspace) COMPLEX array, dimension
  110.              (N) if SIDE = 'L' (M) if SIDE = 'R'
  111.  
  112. SSSSEEEEEEEE AAAALLLLSSSSOOOO
  113.      INTRO_LAPACK(3S), INTRO_SCSL(3S)
  114.  
  115.      This man page is available only online.
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.